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Abstract 

We consider a simple extension to the previously found gravity solution corre- 
sponding to a boost invariant Bjorken plasma, by allowing components that are 
asymmetric under parity flipping of the spacetime rapidity. Besides the question 
whether this may have a realization in collisions of different species of projectiles, 
such as lead-gold collision, our new time-dependent gravity background can serve 
as a test ground for the recently proposed second order conformal viscous hydrody- 
namics. We find that non-trivial parity-asymmetric effects start to appear at sec- 
ond order in late time expansion, and we map the corresponding energy-momentum 
tensor to the second order conformal hydrodynamics to find certain second order 
transport coefficients. Our results are in agreement with the previous results in 
literature, giving one more corroborative evidence for the validity of the framework. 
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1 Introduction and motivation 



Since the advent of the RHIC experiment, dynamics of finite temperature QCD plasma 
has attracted a lot of study in the past few years. One of the interesting features is 
that the plasma produced at RHIC seems to be strongly interacting, and perturbative 
QCD study can't explain certain phenomena observed, such as jet-quenching and the 
small viscosity-to-entropy ratio, etc. Given the situation, gauge/gravity correspondence 
may be of some help to understand at least certain aspects which are common to a wide 
class of examples of strongly interacting gauge theories at finite temperature pLj • A slightly 
simpler version of conformal field theories, the AdS / CFT correspondence, always contains 
a 5D Einstein gravity with cosmological constant in asymptotic AdS^ spacetime, and the 
predictions from it will be universal ones. The hope is that real QCD plasma shares some 
of these predictions. Besides that, it is also of purely theoretical interest to study finite 
temperature plasma dynamics in the dual gravity picture. In the gravity dual picture, a 
finite temperature plasma in equilibrium is described by a black-brane, a still mysterious 
object, whose Hawking temperature is identified as the temperature of the field theory 
plasma [2] , and the small near-equilibrium gravity dynamics on this black-brane spacetime 
reproduces the expected hydrodynamic description of a finite temperature gauge theory 
plasma [3]. The rewards are various transport coefficients of strongly interacting gauge 
theories near equilibrium, which cannot be computed by current field theory techniques 

Although the RHIC plasma is believed to be quickly local-thermalized to facilitate 
a hydrodynamic description [6j, it is not a near-equilibrium plasma; it is an expanding 
plasma with continual decrease of temperature, whose ultimate state would be a cold 
gas of hadrons. Only the late stage of asymptotically slow variation of temperature can 
be described by hydrodynamics. Therefore, the black-brane for a plasma at equilibrium 
is not a priori correct object one can start with, and one needs to look for a full non- 
linear gravity solution to describe an expanding plasma. The early time dynamics of the 
collision, such as isotropization [7] , would require a hard-core numeric analysis of Einstein 
equation [HI El [IDl [Hj, but if one is interested in only late time asymptotic regime, a 
systematic perturbative approach in terms of suitably chosen expansion parameters may 
be invoked to find a solution. This was first initiated by Janik-Peschanski [12] with several 
notable developments afterward in Ref. [El [H Hg [IH [HI [HI [19] . 

A useful assumption in studying the relevant expanding plasma is the boost-invariance 
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first proposed by Bjorken |20j, which simplifies the gravity analysis drastically. It has a 
well-accepted physics motivation from what would happen in the collision of two heavy 
nuclei. We refer the readers to many nice reviews on that, for example Ref.[21j. To 
explain kinematics in few words, take the coordinate transformation to 

= Tcoshy , = rsinhy , (1-1) 

where {x^,x^) are the Minkowski time and the collision direction respectively. The flat 
metric then looks as 

2 

ds^ = -dr^ + T^dy^ + ■ (1.2) 

1=1 

A boost transformation along x^ is simply a constant shift of rapidity y, and the boost- 
invariance implies the homogeneity along y. For simplicity one also assumes independence 
of transverse directions x^'^, which is ok at least in the region nearby central collision 
axis, so that all physical quantities would depend on the proper-time r only. For the 
hydrodynamic plasma, this means that the energy- momentum tensor is written as T^'^(r). 
One more thing, which is again well-motivated in collisions of two identical projectiles, is 
the symmetry under the flipping of spacetime rapidity y, 

P : y^-y , (1.3) 

which is also a space parity transformation x^ —x^. This assumption would put the 
ry-component of the energy- momentum tensor zero; T'^'^^t) = 0, or equivalently the 
local fluid component has a 4- velocity tt^ with u'^ = 1 being the only non-zero component. 
Given these inputs, the corresponding gravity dual solution for the boost invariant plasma 
has been analyzed extensively. 

Our aim is to add a slight twist to the above by allowing components that are odd 
under the parity flipping of rapidity y, while still taking boost-invariance for granted; 
simply put, our plasma of interest will be homogeneous with non-zero T'^^(t) ^ 0. This 
would complicate the local fluid velocity tt^ with a non-zero u'^ 0, whose value is not 
clear a priori and should be determined along the way in the subsequent analysis. We 
stress that the existence of T^^(r) 7^ has nothing to do with the total momentum along 
^3 ^p3^ original Minkowski coordinate which can always be made zero by suitable 

boost transformations going to the center of mass frame. An easiest way of seeing this 
is to recall that boost-transformation is a simple shift in y, under which T'^'^It) 7^ is 
invariant. The P^ is in fact given by 

P^ = f d{Tsmhy) T{l + 2smh'^y)T^y , (1.4) 
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Figure 1: A schematic picture for parity asymmetric boost invariant plasma. 

so that the shift in y indeed affects it, but the statement that there is a region with a boost- 
invariant T'^'^{t) 7^ is something which is not related to whether one is at the center 
of mass frame or not. Rather, the region with T'^^(r) ^ has an intrinsic asymmetry 
that presumably arises from two different species of projectiles, such as lead-gold collision 
for example. See Figure 1 for a schematic explanation. It would be interesting to think 
further whether this is a realistic possibility in experiments l]. 

Besides the question of reality of our parity-asymmetric boost invariant plasma, there 
is another motivation to study it, perhaps more important in the context of AdS/CFT 
correspondence. The gravity solution we will construct is one new time-dependent gravity 
background, whose late time asymptotic behaviors can provide us consistent cross-checks 
for the recently proposed second order conformal viscous hydrodynamics in Ref . [25l [26] . 
Especially, the second order transport coefficients have been obtained by several other 
methods, and it would be a non-trivial corroboration of the framework if one finds agree- 
ment between them and our results, which is indeed the case as we will see. Finally, we 
should mention that our gravity analysis is guided by previous works in a large extent, 
especially those in Ref. [171 [IE] . In most cases, we will follow the conventions and notations 



^ Near the completion of this work, we became aware of Ref.[^ which considered R-charged plasma 
with leading order parity asymmetric term, Ref. [231 which considered a similar construction in 1+1 
dimensions, and Ref.[24J which studied asymmetric collision of two shock waves in AdSs. 



in Ref.[l8]. 
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2 Precursory hydrodynamics 



We will be back to the full second order conformal hydrodynamics in section 5, but before 
presenting the gravity analysis in the next section, it is helpful to consider a preliminary 
hydrodynamics to see what one would expect from the additional parity-asymmetric term 
we are putting. The important thing we will find is that its first non-trivial effects appear 
in the usual parity-symmetric terms at second order in late time expansion. 

The conservation equations D^T'^'^ = in our Rindler-like coordinate (11.21) are 

dr (tT^^) + r'^T^y = , 
dr {tT^^) + 2T^y = , 

drUT^') = , (2.5) 



and the traceless condition is 

- T^^ + T^T^y + T^^ + T^^ = . (2.6) 

One easily integrates the second equation to get the parity-odd term 

4 C 

T-' = --- , (2.7) 

with an integration constant C which we will keep non-zero. We do abandon T"^* however 
in this work. The remaining equations are as usual as in the parity-symmetric cases, 
which one can solve in terms of a single function 

T"" = -a(^)(^) , 

= -i(a(4)(r) + ra.a(^)(r)) , 

= -^(2a^^\T) + Tdra''^\T)) ,z = l,2 , (2.8) 

where the notation is chosen for later convenience. The function a^^^r) contains all the 
non-trivial dynamical details of the expanding plasma, and depends on the microscopic 
theory as well. Note that the parity-asymmetric component T"^^ ^ has a fixed, trivial 
structure ~ due to kinematics, but its effects on a^^^r) may well be highly non-trivial; 
in fact there is no a priori expectation one can make except that it should be an even 
function on C. There is no further information one can extract from the above kinematic 
constraints. 
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One might be misled to think that a^^^r) would be simply decoupled from T^^, but 
it is not true in general as the plasma is interacting non-linearly. To see this clearly, let's 
try to map the energy-momentum tensor into an ideal hydrodynamics form, which is only 
a crude leading order approximation in the late time expansion but would be sufficient to 
make the point here, 

T'"' = e{T)u^u'' + p{t)A^"' + (higher derivative terms) , A"" = gf"" + , (2.9) 

where p(r) = |e(T) for conformal plasmas. For parity-symmetric boost invariant plasma 
with u'^ = 1 being the only component in -u^, one can insert the above into (12.81) to get 
the equation 

e = -^Tdre + --- , (2.10) 
with a solution as the leading order behavior in large r, (we have put Cq = 1 for simplicity) 

e(T)~i, + --- . (2.11) 

T3 

Then consider our situation of having a non-zero T"^^, and an inspection requires us to 
have a non-zero from 

T^y = ^e{T)u^uy + ■ ■ ■ . (2.12) 
3 

Assuming that the leading large r behaviors of e(r) and u'^ are unaffected by the parity- 
asymmetric term, 

e{T)r^X + --- , u^ = l + ... , (2.13) 

r3 

which indeed will turn out to be true consistently in later analysis, one obtains from (I2.12p 
that 

uy^-^ + --- , (2.14) 

ra 

as the first non- vanishing term in the late time expansion. From the normalization — 1 = 
QfivU^u^ = —{u^Y + r'^iu^y, one then finds the sub-leading correction to -u^ to be, 

n" = l + — + ■•• . (2.15) 

2t3 

Back to (12.91) with this u^^, one finds that sub-leading corrections coming from the parity- 
asymmetric term to the parity-symmetric components of the energy momentum tensors 
are ~ -\ smaller than the leading order terms. This is a second order correction in 
late time expansion with -\. Although we need to consider higher order viscous terms 
that we neglected in (12.91) to do a proper analysis at this second order, which will be 
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the subject in section 5, this prehminary analysis shows that non-trivial effects from the 
parity-asymmetric component T'^^ do appear in a^^^r) starting at second order in the 
late time expansion of This input will be helpful to the gravity analysis in the next 
section. 



3 Gravity solution 

Any AdS / CFT set-up allows a consistent truncation of the 5D bulk dual theory to a pure 
Einstein gravity with negative cosmological constant. As there is no global charge in the 
plasma we are considering, we work within this gravity sector only. See Ref . [271 [2H1 [291 [2Q] 
for R-charged plasmas. We take the standard convention of putting the AdS^ radius i = 1 
and the cosmological constant A = —6. This effectively fixes the 5D Planck constant to 
be a specific value depending on the model. For A/" = 4 SYM, it is 

87rG5 = — . (3.16) 

c 

The equation of motion we have to solve for the dual geometry of expanding plasma i 

Emn = Rmn — -Rqain — QQmn = Rmn + ^9mn = . (3-17) 



Our subsequent analysis is largely based on the previous works starting from Ref. |12j 
where the late time expansion scheme of \ was first proposed in the Feffer man- Graham 
Ansatz. It was later found in Ref. [T71 ITS] that an Eddington-Finkelstein form of the metric 
is more suitable for all order regularity of the solution, and this will be our starting Ansatz 
too. Our notations and methods closely follow those in Ref. [18j. 

We propose the following metric solving the 5D einstein equation, 

2 

ds^ = -r^a{r, T)dT^ + 2dTdr + e2(''(-.-)-<=('-'-)) (1 + rrf dy^ + r^e'^''^^ J2 

i=l 

+ 2r f F(r, r) + (1 + ^'^) drdy + - (1 + rr) /i(r, T)dydr , (3.18) 

y 3rr J r 

where the second line is what we add as new terms violating parity symmetry y —y. 
The radial coordinate r is the 5'th holographic direction in addition to the previous 4D 
Rindler-like coordinates (r, y,xi,X2), and r — > cxo is the AdS^ boundary. The factors of 



^Capital letters are for 5D coordinates, while Greek letters are used for 4D indices. 
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(1 + rr) are just for convenience and not essential in late time expansion, but we choose 
to put them according to Ref.jlB]. Near r oo boundary, the induced 4D metric should 
converge to (— dr^ + r^dy"^ + dxf) in order to be asymptotic AdS^, so that the functions 
that appear in the above have boundary conditions 

(a, 6, c, F) ^ (1,0,0,0) , r^oo . (3.19) 

A crucial element in the consistent late time expansion for the gravity solution is to 
introduce a scaling variable 

u = rT^ , (3.20) 

that is kept finite while taking a late time limit r — > oo, and to invoke a power series 
expansion of for every functions to be determined in the metric by solving the einstein 
equation order by order [12]. To be more concrete. 



a(r,r) = ao(M) + + + " " 

7/ N 7 / \ bi(u) b2(u) 

, , , . cAu) C2(u) 

T3 T'i 

F{r,r) = Fo(«) + ^^^") ' 



2 
T3 



h{r,T) = hoiu) + ^ + --- , (3.21) 

where we have shown terms only up to the point of our interest for the purposes of this 
work, but in principle one can go further as he/she wants. It is conceptually straightfor- 
ward to put the expansion (13.211) into the einstein equation (13.171) and to solve order by 
order iteratively. We will find that this works consistently with our metric Ansatz. 

Although it is conceptually simple, it is technically demanding to compute the relevant 
einstein tensors with our metric, and we have used a freeware Mathematica package, 
RIEMANNIAN GEOMETRY & TENSOR CALCULUS [31], for algebraic computations. 
However, what will be important for us at the end is a closed form of the solution up to 
the desired order that we will present in this section. One subtle point which is worth of 
mentioning is that solving Emn = at a given order is not equivalent to solving E^^ = 
at the same order, because various metric coefficients contain different factors of r, so 
that in going from Emn = to Ef^ = 0, the terms at different orders in the former can 
mix up at a same order in the latter. Because the latter seems to be a more invariant 
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notion under the late time scaling r ^ oo (recall that = 5fj ~ for a diagonal 
metric), we will work with it. 

Let us denote the coefficient of t~° in the einstein equation E^' = as E^^^°'\ For the 
zero'th order functions (ao('u), &o(w), Co('u), -Fo(ti), /io(w)), one obtains several non-linear 
equations of them from the leading order terms of the einstein equation. There is no 
definite way of solving these; it is essentially same to finding a solution of non-linear 
einstein gravity. For our purpose, one can proceed by the observation in the previous 
section that the parity-asymmetric terms start to affect parity-symmetric components 
at second order in late time expansion. This suggests that one can simply take the 
previously known solution for (ao(w), foo(w), Co{u)) without parity-asymmetry terms, while 
putting (Fo(u), /lo(ti)) = (0,0) as a leading trial solution; 

,4 



W 



aQ{u) = l-— , 6o(«) = Co(n) = Fo(n) = /io(ti) = . (3.22) 



One easily checks that this solves the leading order einstein equations up to E^^^\ and 
it can be a consistent starting point to do subsequent series expansion. Once the zero'th 
order solution is given, the problem becomes much more easier and systematic, because 
the next leading einstein equations give linear differential equations for the next order 
terms in the solution {ai{u), hiiu), ci(-u), -Fi(m), hi{u)), and so on. One can in principle go 
on further order by order iteratively. 

One finds that the next order einstein equations, Ep'^\ Eu^^\ Et'^^\ Eu^^\ Ey'^'-^\ and 
E^ ^ provide complete linear differential equations for (ai(u), f)i(u), ci(u)) only, without 
involving (Fi(m), hi{u)) at all. Therefore, their solution is simply identical to the previous 
results of the parity-symmetric case, 

, , 2w^{u + w) + Ci{u^ + w^) 

ai{u) = , 

3 

hi{u) = ^ , (3.23) 



ci{u) = — (-^ + 4(l + Ci)- + 2tan-i f-) +log 
Qw \ u \wj 



where Ci is an integration constanllfl which is actually a pure coordinate reparametrization 
freedom, and we have already determined another integration constant to have a regularity 
dX u = w which is the only dangerous point except the trivial singularity at m = hidden 
behind event horizon. This is indeed in line with our expectation from the previous 



^In comparing with Ref.[18|. Ci = — — 1. 
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section that non-trivial effects of parity-asymmetry would emerge starting only at the 
second order. 

One obtains differential equations for (Fi(w), /ii(-u)) from the parity-odd components 

t(--) y{-) y(-) 

of the einstein equation. From Ey , Eu , and Er ^ we get 

uF[\u) + 5F[{u) = , (3.24) 
whose solution with the boundary condition fl3.19p is 

Fiiu) = -, . (3.25) 

u(-) 

The remaining component, Ey ^ , gives us an equation for hi{u), but only in the combi- 
nation of 

Si{u) = hi{u) + ^u^F^{u) . (3.26) 

In fact, if we look at the next order parity-odd einstein equations, Ey^''\ -E^*-^-*, Et^^\ they 
are equations for Si{u) without separate hi{u) or F2{u). All these equations are solved 
uniquely by 

u {21u + lOCi) {u" + uw + w'')+h{h + 2Ci) 

(m + U)) (m2 + U)2) ■ ^ ' ' 

The reason for this combination can be traced back to the fact that the F2{u) can be 
simply gauged away by the coordinate transform 

y-,y+^F2iu)\ + --- , (3.28) 
4 7-3 

under which hi{u) becomes hi{u) + ^u'^F2{u) = Si{u), so that only Si{u) is a gauge 
invariant combination at this order. This means that one can choose the gauge F2{u) = 0, 
which gives finally 

It seems to us that we can always gauge away Fn{u), n > 2 with similar coordinate 
reparametrizations of y order by order, so that there is a gauge with F{u, r) being 

F{u,t) = ^\ , (3.30) 

U T 3 

exactly without any higher-order modifications. This seems to be consistent with the fact 
in the previous section that T'^^ has a fixed structure of without modifications due to 
kinematics. 



9 



Having obtained {Fi{u), hi{u)), wc then consider the next order einstein equations for 
{a2{u), b2{u), C2{u)), and one indeed discovers non-trivial effects from {Fi{u), hi{u)) at this 
order. Et ^ , E^^^\ Et ^ , Eu ^ , Ey ^ , and E^ * provide complete differential equations 
for (02(^)1 ^2(^)1 C2(^x)) which include "sources" from various quadratic combinations of 
{Fi{u) , hi{u)) , in addition to the usual terms from parity-symmetric components. One 
can also check that there are no further contributions to these equations from higher order 
terms in the expansion. We simply present the solution of these equations after fixing 
certain integration constants to be regular at u = w, 



a2{u) 



+ 1 + ^ (6 + Ci) - 2C2«^ -il + C,f—-2w'(l + ^ + C2w) 
6 u \ 6 J 

iu' + w") tan-^ ( -\ + u^^ + -u'W + -u^'^w^ + -^i^w' + -uW 



^ + w') tan-^ ( -) + v}' + \u'W + lu^'w' + ^\'w' + \ 

„i6^6 J \yjj 3 5 45 9 

,9 



~ 308 ^ It + " ^ ~ 7^ J \ [u + w)\ , (3.31) 



' ^+^r-i+c,ri+^)+3C2« 



3uw \ u \ V 2 
^2/ 3 /-u\ ^ ^ /_3 ^ ^ 3 ^ \ 

75 (5 + 2C^f ( luA \ 3 (5 + 2C,) \ 

224 u8 llu')^ 8 u^^) ' ^^-^^^ 



c'Ju) = ^— J- 1 (u^ + 2mw^ - 3i(;^) tan"^ f-) + u(u^ log 



-\- log (^^ ~^^Mm ~ ( TTU^ — U^W — UpW^ 

- (3 + 7r)u^w^ - (7 + 27r + 10Ci)u'^w'^ + 2(2 + Ci)u^w^ + 2(2 + Ci)u^w^ 

+ (2 + 2Ci + 7t)uw'' + (2 + 4Ci + Tr)w^ + (tt - 2C2 - 6C2m)«j(w - w) (^m^ + 2'u^'u; 

+ C'(^ tan-i ( -] + -r:--A ^ ( u'' + 2uV + lu'w' - ^u'w' 



^u^w'^ \wj u^^w^{u^ + w"^) \ 5 70 

5485 4 „ 17315 2 m 50^ , ^ , « / « 42 21 2 A 
~m ^ "396~ ~ y ^ V ^ ~ 22 j 
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^(5 + 2C^)uw\u' + + 1^(5 + 2C^)V ) ) , (3.33) 



where C2 is a constant of integratioiu, which is a pure gauge. Note that the pieces 
proportional to are the new terms coming from our parity-asymmetric components. 

Up to this order, we find the unique regular solution except Ci and C2 which corre- 
spond to a coordinate transformation as observed in Ref.[18]. 

u ^ u + — , (3.34) 

3r3 3t3 

and there is no further ambiguity left in the solution. It seems very plausible that all 
order regularity in the present expansion scheme can also be proved. 



4 Holographic renormalization and energy momen- 
tum tensor 

Our next obvious step is to find the 4D energy-momentum tensor of the parity-asymmetric 
expanding plasma corresponding to the gravity solution we constructed in the previous 
section via AdS/CFT correspondence. This involves the well-established procedure of 
holographic renormalization [321 [S3], and we will be very brief in explaining it. Consider 
a constant, large r hyper-surface as a boundary of our asymptotic AdS^ spacetime, whose 
induced metric we denote as 7^,^. Let K^^ be the extrinsic curvature of the hyper-surface. 
The "bare" energy-momentum tensor computed from the well-known Brown- York method 
for the spacetime with boundary, 

Tlr = ^r\K,^-K^,^) , (4.35) 
turns out to be divergent in r 00 limit. The factor in front is from the fact that 



It^u = r^g'ilu^, so 



= r4-=-'^7^ • (4-36) 



This divergence has a field theory interpretation of UV divergences one encounters for 
the "bare" operators before regularization and renormalization. Therefore the question 
is what are the correct counter-terms one should subtract from the above to have a 
well-defined renormalized expectation value of T^j^. The important ingredient in choosing 



^It is similar to ^2 in Ref. 
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these counter-terms is the requirement of general covariance to correctly reproduce certain 
conformal anomalies. The upshot is that the counter-terms should be local covariant 
objects constructed purely out of the induced metric 7^1, on the boundary. In a "minimal 
subtraction scheme", one choose them minimally just to be sufficient to remove all the 
divergences in the bare terms. The result is 

= ^^lim^^oo?'^ (k^^ - K-f^^ - 37^^ + l;G^J) , (4.37) 
ottGs V / / 

where G^i, = Rfj,u — ^Rj^u is the einstein tensor of 7^1^. 

It might look rather straightforward to implement this procedure by simply plugging 
the gravity solution of the previous section into 04.371) . but in reality one finds it to be 
impossibly complicated to do in this way. A better strategy that has been used in practice 
all the time, and also in the formulation of holographic renormalization itself, is to first 
solve the 5D einstein equation near boundary r ^ 00 in general, reducing the number 
of independent degrees of freedom, and to use these results in computing the energy- 
momentum tensor (14.371) . One then simply needs to identify the relevant near-boundary 
terms, that appear in the energy-momentum tensor expression, from the actual gravity 
solution to get final results. In obtaining near-boundary form of the metric that solves 
the einstein equation, we will use one fact from the previous section; note that Fi and hi 
in our solution have ^ ^ behavior near r ^ 00, so that up to our desired order F and 
h can be safely assumed to have expansion starting 

F(r,r) = — + , h{r,T) = —y + ... . (4.38) 

We also checked independently that these are consistent with the einstein equation near 
boundary. Because the parity-symmetric components a(r, r), 6(r, r), c(r, r) are even func- 
tions on {F,h), the effects from {F,h) to {a,b,c) can be at most ^. As the energy- 
momentum tensor is sensitive only up to ^ terms, these effects won't affect the near- 
boundary structure of (a, b, c) that are relevant for the energy-momentum tensor. We 
stress that this does not mean that the values of (a, b, c) up to ^ near boundary, and 
hence the energy-momentum tensor result, is not affected by the presence of [F, h). Only 
the near-boundary structure is insensitive, whose meaning will be clearer in a moment. 
By assuming (14.381) . explicit computations give us 



a(r,r) = 1 + ^ + 

r 



r 



V 
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^2{b{r,T)-c{r,T)) 



a«(r) 

1 + — + 



( (a«(T) 



+ 



r 



2r 



+ 



3c(r,T) 



(aW(r))^ 



+ 



(4.39) 



a«(r) 
1 + — + 



-V(r)-^r9.a(^)(r))l + .. 



which is essentially identical to those in Ref. [T8] without {F,h) as mentioned before. 
However, the function which is undetermined by boundary analysis does depend 

on (F, h) or equivalently C. Finally, h^^^r) can be given in terms of F^^\t) by einstein 
equation, but we won't need it for our purposes. 

With the above, the renormalized energy-momentum tensor is easily found to be 



T 

J. q-- 



ry 



T 

J. 7 



yy 



T- 

7,7, 



SttGs 
1 

1 

SttGs 
1 



3 



(aW(T) + ray^)(r))) 



4 



1,2 



(4.40) 



Note that doesn't appear because it can be removed by a coordinate transforma- 

tion. From the explicit solution in the previous section, one can easily obtain the required 

a^'^^r) and F*^^^(r) by recovering r = Mr~i, 

;i + 21og2)u;6 + 12C2 
F^'\t) = % , (4.41) 



where we can see the parity-asymmetric effects of C at the second order expansion. 



5 Mapping to second order conformal viscous hydro- 
dynamics 

In this section, we will try to map the energy-momentum tensor from the gravity side in the 
previous section to the recently proposed second-order conformal viscous hydrodynamics 
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in Ref . [251 126] . We remind the readers of that this is not a priori guaranteed to work 
because the proposed hydrodynamics framework is restrictive in describing the energy- 
momentum of plasma, with certain finite number of transport coefficients. Any success 
in this mapping will be a corroboration of the validity of the framework, and we think 
this is at least worth of checking explicitly. We will find that the mapping indeed works 
with some second order transport coefficients identified along the way, whose results are 
in perfect agreement with those in the literature via other methods/backgrounds. 

Simply put, the second order conformal viscous hydrodynamics (SCVH) in Ref. |25] is 
a statement about the energy-momentum tensor of plasma, sometimes called constitutive 
relations, expressed in terms of a finite set of local thermodynamics quantities. In our 
case at hand, these will be 4- velocity tt^ with m^m^ = — 1 and the energy density e. The 
pressure p = |e is dictated by e by conformal nature. There are four unknowns and we 
have four equations from D^T'^'^ = 0, so that one is given a complete, closed system of 
dynamics. The proposal was 



(Da""^ + ^a"" (V ■ u) 



+ K 



i?^^'^^ - 2uc,R"^>"'^'^Uf^ + Xia^^a^'^^ + X2a\n'''^^ + X^n^^n'^'^^ , (5.42) 



where A^'^ = g'^^ + u^u^ is the projection to the transverse space to u^, D = u^V and 

^m- = 2V<'^m"> , n^"" = ]^A^''A''^ {Va.up-Vpu^) , (5.43) 

with 

^ 1 ^M^A-/^ (A„^ + Ap^) - ^A^^-A^^^A^p , (5.44) 

for any two-tensor A^^ . The first order coefficient rj is the shear viscosity, while r//, and 
-^1,2,3 are second order transport coefficients. As our 4D metric is fiat and the vorticity 
Vt^^ vanishes for our boost invariant plasma even with parity-asymmetry, we won't have 
an access to k and A2,3 by our present expanding plasma. 

The quest is to choose the right e(r) and u^ij) as well as r], tu, and Ai to reproduce 
the energy-momentum tensor from the gravity side we obtained in the previous section. 
Since the gravity result already solves D^T^'^ = 0, this automatically includes the conser- 
vation equation in the hydrodynamics side too. In fact, we would need one more piece of 
information; a general expectation is that local thermodynamics is completely specified 
by e onljl^l, so that the transport coefficients should also be determined completely in 



^or equivalently by the local temperature T which is related as e = ^^N'^T^ for A/" = 4 SYM 
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terms of e only. By dimensional counting, one has in general conformal plasma, 

rj = r]^e^ , Ai = A?e5 , m = r^je-^ , (5.45) 



with fixed dimensionless numbers r]^, A^, and rfj, which are intrinsic to the microscopic 
details of the theorjj^. At the end, the aim becomes : choose e(r) and tt^(r) as well as 
fundamental dimensionless constants 77°, A^, rfj to match the energy-momentum tensor 
from the gravity result. The non-trivial test of the formalism would be that r]^, A^, and 
Tjj should be same to those obtained from other methods/backgrounds. 

From the discussion in section 2, we know that e(T) and 'U^(r) start their expansion 

as 

e(^) = TTT^ — + — + — + ■ 



T3 I rs 



u'{r) = — + 4 + ...)+ O(C^) , (5.46) 



where we restrict ourselves up to 0{C^) because our results in the previous section f l4.4ip 
can tell things only up to this order. The u'^ is given from w^m'^ = — 1 as 



U T 



V^lT72|i^K^ • (5.47) 



It is computationally straightforward to insert the above general expansions into fl5.42p to 
have an energy-momentum tensor from SCVH, and to compare it with the gravity result 
(14. 40 p . (14.411) order by order to see whether it works. 

From comparing T^^, one has 



and from T^y, 



3wei + 2 = {t-') , (5.48) 
(l8z/;2e2 - 1 - 2 log 2) + 12^2 = (r"t) , (5.49) 



/167rG'5\3 r/° 



2/1 + ei - (^^-^ j ^ = (r-3) , (5.50) 



where one can check that the next order r 3 involves 0{C^) effects for which we can't 
say much. From the comparison of Tyy, we have 



1 



2 + u;e,-4(^^^)%o = (rO) , (5.51) 



^These definitions are different from those in Ref. I 
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18w\2 - (1 + 2 log 2) - 54 — ^ ) t/j'^eir]' 



2 , 



+48(1:::-^) ^2 AO_^o^o^^— ^ Q ^^-t) ^ (5 52) 



and finally from Tu [i = 1,2), 

«;ei + 2^-^ = (r-2) , (5.53) 



18^^62 + u;^ (1 + 2 log 2) + 27 (^^) ' wh,r]^ 

-24(l^)^^(A?-A{^)+i^ = (r-f) . (5.54) 

One has five unknowns ei,2, Hi, if, (A° — rpT^i) with the above seven equations to solve, 
and it is not a trivial thing for the SCVH to work, but one indeed finds the consistent 
solution to the above, 

2 _ (l + 21og2)w;'5-12C2 _1 

To see the above results are in agreement with the literature, it is convenient to rewrite 
things in terms of local temperature T by using 

so that (I5.45P and ( I5.55P give us 

' 16^5 ' ^ 1 ' 2 I6G5 

Recalling that the area element at the horizon of a static black-hole with temperature T 
is 

r% = -K^T^ , (5.58) 
so that the entropy density is related to the shear viscosity as 

which is a famous result. The Ai and th that are obtained in Ref. [25] for A/" = 4 SYM 
are 

A, = JL , ,„=(i^i^ , (5.60) 

which are again in agreement with fl5.57p . Our main point is that the mapping works fine 
even after our parity-asymmetric components included. 
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6 Discussion 



We present a simple extension of previously studied boost invariant plasma by allowing 
parity asymmetric components, and obtain its corresponding late time gravity solution 
up to second order expansion in the context of AdS/CFT correspondence. We check that 
the second order conformal viscous hydrodynamics in Ref.[25] consistently describes this 
plasma. There are a few directions that may be worth of pursuing further; we haven't 
looked at the location of apparent horizon and the entropy density to see whether the 
claim in Ref. [34] is valid in our new background. Another direction would be to relax 
assumptions in the gravity solution further to find more general gravity solutions, for 
example including radial profile and/or anisotropicality in the transverse plan^. Numeric 
analysis for early time dynamics within these more general setting in AdS / CFT may also 
be pursued further. 
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A Appendix 

In the appendix, we list the necessary expansion of the einstein equation up to second 
order, after fixing the zero'th order solution 

4 

W 

ao{u) = 1 J , bo{u) = co{u) = Fo{u) = ho{u) = . (A. 61) 

The prime denotes derivative with respect to u. 



El 



2^5^2/3 

+ Au'hi - 2Au^F'l - ASu^Fihi - 2Au^hl + I2u^w%l - I2w^h\ + ^u^a[ - 2u%[ 
- Qu^aib[ - 3u^^a[b[ + Qu^c[ - 21u^F^F[ - 6u^h^F[ - Qu^w^hF[ - 3u^^F[^ 



^See Ref . [351 ESI EZl |38j for some studies on anisotropic situations, and Ref.p9] for including radial 
profile. 
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+ 
+ 



—2h' — vh" —2h' — vh" 1 / 
" mtV3 ut 2uH\ 

- 3w^cf - u^h^F[ - 6u^Fih[ - 6u^hih[ + 6uw'^hh[ + u^hF^^ , 

+ 2u^b'. - 2u^w% + 3?/c' - Zv:'w^cA + -^-—Il2u^w^ - Su^ai - Au%i 

+ 6u^w\ - 6u^ci - 6u^w% + ISu^F^ + 36u^Fihi - 36uWFi/ii + ISu^hl 

- Seu'^w'^hl + ISw^hl - 3u^a[ - 3u%ia[ + 8u%[ - 6u'^w%[ + 3u^aib[ 
6u\b[ - 6u^w%ib[ - 6u^cib[ + 6u^w%b[ + 6u^w% + 9u^aic[ - 6u%ic[ 

El = 2Jt^/^ (^^^ ~ ^"^"^ ~ ^^^"'^ ~ ^"^^'^ ^ '2u^w%\ - u'^a'l - 2vJ\J[ + 2v?w'^\l^) 
+ 1 .1. f-Sw^aa - 7M^a' - ^u^% + 2u^w% - u^ol' - 2u%'' + 2u'^wX) 

+ } A,'. I - Qu^w^ + 6M^ai + Wbi - 24u^F^ - AM^F^hi - 2Au^h\ 

2Au^w^hl + 2Aw^hl + 3u^a[ + 2u%\ - 12uWb[ - 18u\ib[ - 3u^'^a[b[ - 6u^%f 
6u^w%'^ - 6u^c[ + 12uW, + 12u^%[c[ - \2u^w%(!^ - 9ii^°cf + 9u''w^(!l 
A2u^FiF[ - 2Au^hiF[ + 12u^w%F[ - 3M^°Ff - 2Au^Fih[ + Uu^w^^F^h'^ 

- 24u%ih[ + 36u^w*hih[ - 12uw^hih[ - 4u%'i - 6u^\ib'l - Gu^^F^F^^ , 

El = — f-4«3ai - u^a\ - 6u%[ + 2w% + 5u^c[ - w^c[ - u% + uw%'[ + U'd! - uw^c 

+ 3./^ ( 24«^ai + 16m^6i - \2vJ - 24u^F^ - 48u^Fihi - 24u^hl + Uu'^w'^hl 



+ 

- 42^9 



+ 
+ 



I2w^h\ + Gu^a\ - Au%\ - I2u^w%'^ - 3Gu\ib'^ - Gu^^a'^b'^ - &u^%'^ + Gv!'w%'^ 
. 2m^c; + e>u^w^c\ + 3Qu^aid^ + Q>u^^a\d^ + Q>u^%\d^ - Q>u^w%\d^ - 27u^FiF[ 

- Qu^hiF[ + 6u^w'^hiF[ - 3ii^°Ff - 6u^Fih[ + Qu^w*Fih[ - Qu^hih[ + 12u^w^hih[ 

- 6uw^hih[ - 4u%'l - Qv}^aib'[ + Au^(J[ + GK^^aic'/ - 3v}^F^f{^ 
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+ \2u^w%\ + 12w^/i? - 2u^\)^ - 6M^ai6; - m^c^ - 'iu^w^d^ - \^u^a^d-y - 'iu^^a^d^ 

- 3m^°6'ic; + 'iu^w%d^ - Uu^FiFi - 6u^hF{ + Qu^w%F{ - Qu^F^h^ + 6u^w^Fih[ 

- Qu^hih[ + 12u^w^hih[ - Quw^hih[ - 2u^d[ - Su^^aic'/^ , 

El = ^[5uF{ + u^F[')t'^'' + -^-^hFi + 12hi-30u^F,b[ + 30u^Fid, + ^^^ 



3u%F[ + 6u^d^F[ + Ihu^F!^ + Auh\ - Qu^Fih'[ + Qu^Fid[ + 3u^F'l + Sm^F. 



2 



-5tx^F{ + 5^/;^F{-i.5F,-- + ^ti;X , 1 ( 194, ^..4, 

= ^^^3^575 + ^^71:^(^-9^ ^1 + 3^ ^i-12u/^i + 12u;/ii 

+ 15u^Fia[ - 12uWFi6; + 24u2«;^Fic'i + 8u^F[ - ISuw'^Fl - Ibu^a^Fl + 30u%iF[ 

- 30u^w%iF[ - 30u^CiF[ + 30u^w%F[ + 3u^h\F[ - 3u^w%\F[ - Qu^d^F[ + Qu^w^d^F[ 

- Ibu^F!^ + Ibu^w^F!^ - Au^h\ + 4Mti;^/i; + 3u^ Fia'i + li^Fi" - 3u^w*F'l - 3u^ a^F'; 

+ e^/^foiFf - ^v?w%rF'( - Qu'ciFl + Qu^w^CiF'l - 3u^F^' + Su^u^^F^' ) , 



/7" _|_ 7; /?" 



- 3u^h\F[ + 6m^c;F{ + Ihu^F!^ + 4m^/i; + Sw^/iia;' + Qu' Fid[ + 6m^/iic'; - Qu^w^hid[ 

- Sti'^F" - 6w^6iF;' + 6w^CiFi" + 3m^F2 ) . 
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